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^-h '. New simulations demonstrate that low-mode, nonradial hydrodynamic instabilities of 
the accretion shock help starting hot-bubble convection in supernovae and thus support 
explosions by the neutrino-heating mechanism. The prevailing conditions depend on 
the high-density equation of state which governs stellar core collapse, core bounce, and 
neutron star formation. Tests of this sensitivity to nuclear physics variations are shown 
for spherically symmetric models. Implications of current explosion models for r-process 
nucleosynthesis are addressed. 

^ ■ 

Or 

I. introduction 

h ; 

Supernova (SN) models in spherical symmetry (ID) — even the latest ones with "Boltz- 
mann neutrino (u) transport" [Q] — do not explode, neither by the prompt bounce-shock 
mechanism, nor by the delayed neutrino-heating mechanism. Suggestions that neutron- 
^ ■ finger convection in the neutron star (NS) could enhance the v luminosity and thus the 
heating behind the shock [|2] could not be confirmed by detailed analysis Also Ledoux 
convection in the nascent NS seems to be ineffective in this respect [EJII3- A reduction of 
the v opacity by several 10% at densities below ~ 10 13 g cm -3 would also allow z/'s to leave 
the accretion layer more easily and could produce explosions [ E] , but a physical effect of 
corresponding size is not known. 

Multi-dimensional simulations of core-collapse SNe revealed that hydrodynamic over- 
turn behind the stalled accretion shock develops on a relevant timescale and increases the 
efficiency of the //-heating mechanism, enabling explosions even when spherical models 
fail [Ej- Initially there was hope that this new twist in the theory of SNe was the long- 
sought guarantee for "robust" explosions. But the first generation of mult i- dimensional 
models employed a simplified treatment of the v transport — at best by grey flux-limited 
diffusion schemes — which fell much behind the sophistication of the transport in non- 
exploding ID models. It was therefore suspected that the transport approximations might 
have contributed to the success of the multi-dimensional simulations [|S|. 

Recently, two-dimensional (2D) models have become available with a significant im- 
provement of the v transport by solving the energy-dependent equations of v number, 
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energy, and momentum, making use of closure relations obtained from a model Boltz- 
mann equation [ |U| . The description of //-matter interactions has thus reached a new 
level of accuracy and refinement, although the transport is still treated to be essentially 
radial (lateral flux components are ignored in the moments equations but terms associated 
with lateral velocities and lateral gradients of v pressure are included). Simulations with 
thus improved treatment of the v physics could not find explosions despite of convective 
activity in the cheating layer behind the SN shock confirming the suspicion that 
radical transport approximations might have favored the rapid and powerful explosions 
seen in other 2D and 3D models [17} HUj. 

So what is missing in the currently most advanced SN models which fail to produce 
explosions? In the following, we shall first briefly summarize the status of core-collapse 
modeling by the Garching group. In particular, we shall discuss that any limitation of 
the angular wedge of 2D simulations to less than 180 degrees (e.g. to 90 degrees as in 
many previous simulations) imposes artificial constraints to the fluid flow and suppresses 
large-scale modes that can play an important role for the growth of convection. It is 
further demonstrated that even a "modest" amount of rotation in the SN core, in the 
ballpark of predictions of current stellar evolution models [ 111), may have an important 
impact on the postbounce evolution. Finally, we shall elaborate on uncertain aspects of 
the nuclear physics which the SN shows sensitivity to. 

2. EXPLOSION MODELS 

Simulations leading to explosions can be reported for stars near the low-mass end of 
SN progenitors, namely for stars in the mass range of ~8-10 M & with O-Ne-Mg cores [ 
IT!?] , and for an 11.2 M & progenitor [Eh both characterized by small cores of less than 
^1.3M (O-Ne-Mg or Fe, respectively) and a steep density decrease (entropy increase) 
outside. 

2.1. Stars in the 8— 10 M Q range with O-Ne-Mg cores 

The main improvement of our new simulations of O-Ne-Mg core collapse — which 
we did so far only in spherical symmetry - - compared to previous approaches is the 
more accurate, spectral treatment of v transport and //-matter interactions. Using the 
nuclear equation of state (EoS) of Lattimer & Swesty [ HI] and more recently also that of 
Hillebrandt & Wolff [U^j , we could not confirm the prompt explosions found in calculations 
with simpler v treatment [ I15j. The shock is created much deeper inside the collapsing 
core than in the "old" models (cf. also Ref. [HE]), typically at a mass coordinate around 
~0.45 M & [EI (defined by the moment when the postshock entropy first exceeds 3 fee per 
nucleon), and it stalls (defined by the time when the postshock velocity becomes negative) 
only 1.2 ms later at ~0.8 M Q , still well inside the neutrinosphere and before energy losses 
by the prompt v e burst could have contributed to its damping [IT7]. We also do not find 
the powerful shock revival by v heating as seen in Ref. [|ll| and for some choice of the 
nuclear EoS by Fryer et al. [ HH] . Instead, the shock continuously expands due to the 
monotonically decreasing preshock density and the steep density decline at the interface 
between C-0 shell and He shell (Fig. Q). At the end of our simulation the mass accretion 
rate by the shock has correspondingly dropped to less than 0.03 MqS" 1 . Although our 
simulation is not yet finally conclusive in this point, we see indications that a z/-driven wind 
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Figure 1. Left: Composition in the highly degenerate 1.38 M© core of an 8-10 M progenitor [ 
IT2] with O, Ne and Mg in the central region enclosed by a C-0 shell. (X) denotes the mass 
fraction of a representative neutron-rich heavy nucleus that appears in nuclear statistical equi- 
librium, and Y e is the electron fraction. Right, top: Mass trajectories, shock position (thick solid 
line that rises to the upper right corner of the plot, where it reaches the surface of the C-0 shell), 
neutrinospheres {y e : solid; v e : dash-dotted; f M , u^, v T , v T : dashed), and gain radius (dotted) for 
the collapsing O-Ne-Mg core. The mass trajectories are equidistantly spaced with intervals of 
0.02 Mq. The outermost line corresponds to an enclosed mass of 1.36 Mq near the outer edge 
of the carbon shell. The (blue) hatched region is characterized by a dominant mass fraction 
of alpha particles. Right, bottom: Luminosities of u e (solid), v e (dashed) and heavy-lepton z/'s 
(individually; dotted) as functions of time. The sawtooth pattern during the collapse phase until 
about 50ms is a numerical artifact. (The plots were taken from Ref. [I17j.) 



begins to fill the volume between SN surface and shock and will lead to mass ejection with 
a rather low explosion energy (a few 10 50 erg due to the wind power and nuclear burning). 
Little nickel production (~0.01M o ), can be expected, corresponding to the wind mass 
loss rate that can be estimated for an initial v e + p e luminosity of ~ 8 x 10 52 ergs -1 (Fig.Q) 
using the equations in Refs. [1201- The baryonic mass of the NS will be very close to or 
only some 0.01 M & less than the C+O core mass (which is 1.38 Mq). These findings are 
very similar to the outcome of simulations of accretion induced white dwarf collapse to 
NSs (AICs) [EH- 

A weak explosion of an 8-10 M Q star has been suggested as an explanation of the 
observed properties of the Crab SN remnant [1221- The wind-driven explosion seen in our 
models, however, does not provide the conditions for the low-entropy, low-F e r-process 
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Figure 2. Three stages (at postbounce times of 141.1ms, 175.2 ms, and 225.7ms from left 
to right) during the evolution of a (non-rotating) 11.2 progenitor model from Ref. [ I13j. 
visualized in terms of the entropy. The scale is in km. The dense NS is visible as low-entropy 
circle at the center. The computation was performed in spherical coordinates, assuming axial 
symmetry, and employing the "ray-by-ray plus" variable Eddington factor technique [ |2U f° r 
treating v transport in multi-dimensional SN simulations. Equatorial symmetry is broken on 
large scales soon after bounce, and low-mode convection begins to dominate the flow between 
the NS and the strongly deformed SN shock. The model continues to develop a probably weak 
explosion, the energy of which was not determined before the simulation had to be stopped 
because of CPU time limitations. 



nucleosynthesis discussed for prompt explosions of collapsing O-Ne-Mg cores in previous 
work [1231- R- processing in the high-entropy environment of the ^-driven baryonic wind 
can also not be expected to take place, because sufficiently high entropies, short expansion 
timescales, and low proton-to-baryon ratios require the NS to be very massive (~2 M Q ) 
and compact (^10 km) (e.g., Refs. [ 20 ). It is therefore unclear how O-Ne-Mg core 
collapse events could contribute to the production of high-mass r-process nuclei. 

2.2. Massive stars with iron core 

Progenitor stars in the mass range between 11 M Q and 25 M & were found to neither 
explode by the prompt bounce-shock mechanism nor by the delayed ^-heating mechanism 
in spherically symmetric simulations [ I27j . in agreement with results of other groups [Q]. 
This implies that state-of-the-art SN models do not support suggestions that r-process 
nucleosynthesis might occur in prompt explosions of low-mass (~11 M ) progenitors [ 28J. 

We have also started to perform core-collapse simulations for the mentioned stellar 
mass range in 2D, using a polar coordinate grid and assuming azimuthal symmetry [HUE]- 
Because of the energy-dependent treatment of v transport and z/-matter interactions the 
requirements of CPU time are substantial, and we were so far able to perform only a 
handful of such 2D calculations. To limit the need of computer resources we initially 
constrained the computational volume to a ~90° wedge (between roughly +45° and —45° 
around the equatorial plane and periodic boundary conditions for nonrotating models 
and between 0° and 90° with reflecting boundaries for rotating ones), using an angular 
resolution of ~1.4°. This choice was also motivated by 2D models of the first generation 
which were able to obtain explosions with a similar setup due to the help of hot-bubble 
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Figure 3. Maximum and minimum radii 
of the shock (giving a rough measure of 
deformation) vs. time for 2D simulations 
with a polar coordinate grid (models "2D"), 
compared to ID results with the same 
physics, but no convection (models "ID"). 
Top: 11.2 Mq star [ ^], simulated with 
a full 180° grid (Model 2D(128)) and with 
a 90° angular wedge around the equatorial 
plane (Model 2D(64)). Bottom: 15 M star 
(Model sl5s7b2 of Ref. [EUD without (Model 
2D(32)) and with rotation (Model 2D(rot); 
see also [El EE])- In the former 2D simulation 
a 90° wedge around the equator was used, in 
the latter a 90° grid from pole to equator. 



convection [E|. 

First simulations of stars of 11.2, 15, and 20 M & with a 90° wedge, however, did not 
produce explosions [E|. This suggests that the success of the previous calculations with 
simplified transport was most probably connected with the use of transport approxima- 
tions. One of our new models, a 15 M star (Model sl5s7b2 of Ref. [EE]), was a ^ so 
collapsed with a modest rotation: The pre-collapse core was assumed to rotate (essen- 
tially rigidly) with a period of 12s, i.e. Q ~ 0.5rads _1 (see Fig. 1 in [EE!), which is 
significantly faster — but not orders of magnitude faster — than predicted by the latest 
stellar evolution models [E]- Even this "modest" amount of rotation turned out to make 
a big difference. As visible in the lower panel of Fig. El the shock expands to a much 
larger radius than in the 2D simulation without rotation (although both simulations did 
not lead to explosions within ~270ms of postbounce evolution), mostly because of the 
influence of centrifugal forces and the more violent convection in a more extended region 
of v heating [EJ- 

2.3. Nonradial shock instabilities and low-mode postshock convection 

While the 11.2 M & simulation with 90° wedge did not explode (see Fig. El upper panel), 
the same model with a full 180° grid developed a presumably weak explosion (Figs. Eland 
EJ). Equatorial symmetry is broken on large scales some 10 ms after convective activity 
in the postshock region has started (at ~50ms p.b.), and low modes (I = 1, 2) begin to 
dominate the flow pattern between NS and shock after about 140ms p.b., a phenomenon 
which might be linked to the observed large recoil velocities of young pulsars [EH]- The 
convection becomes significantly more violent than in the 90° simulation where obviously 
important degrees of freedom were suppressed. This latter fact was emphasized in an 
interesting paper by Blondin et al. [EDI, who found a similar development of low-mode 
asymmetries in the postshock flow and shock oscillations as a consequence of the instability 
of standing accretion shocks ("SASI") against non-radial perturbations that are amplified 
in a "vortical-acoustic feedback cycle" [ 131 j . 

Blondin et al. considered idealized conditions in their numerical studies, assuming 
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steady-state mass accretion with a constant rate, a fixed inner boundary radius, a simple 
ideal-gas EoS, and ignoring v heating and cooling. We decided to test the effects of z/'s 
and of the size of the angular wedge (and thus of the possible modes) in a separate set 
of simulations with a realistic EoS but an approximative treatment of the v transport (as 
described in Refs. [1211), which allowed us to save computer time and thus to perform 
faster calculations with higher resolution, and to run more models. The mass accretion 
rate was given by the collapse of a 15 M Q progenitor star and the NS was replaced by a 
gravitating point mass inside a contracting inner boundary that followed the behavior of 
the shrinking NS in our full-scale SN simulations with detailed transport physics. 

With this setup we confirmed that 2D simulations with a 180° grid can yield explosions 
even if models with 90° wedge do not explode. In these studies, with v effects switched 
on or off, we found that v losses promote the action of the vortical-acoustic cycle by 
allowing matter to settle on the proto-NS surface. Our results show that corresponding 
shock deformation produces growing perturbations in the postshock flow which can accel- 
erate the onset of convective overturn even in cases where the infall timescale is initially 
shorter than the growth timescale of Ledoux convection. We also see that dipolar shock 
oscillations become more violent in case of a rapid contraction of our inner boundary 
(mimicing a NS that becomes rapidly more compact due to a soft nuclear EoS), thus re- 
leasing gravitational binding energy which partly is converted to turbulent kinetic energy 
of the postshock flow. In this case 1 — 1,2 modes begin to dominate faster and the SASI 
conspires with convective instability to establish favorable conditions for a high efficiency 
of v heating behind the shock. Both nonradial instabilities differ characteristically in the 
way how anisotropies emerge: Convective overturn is driven by the negative entropy gra- 
dient in the z^-heated layer, which becomes Rayleigh- Taylor (RT) unstable first on small 
scales, while the SASI starts from vorticity producing sound wave interactions with the 
shock and induces dipolar oscillations of the postshock layer before convective activity is 
initiated and the typical RT mushrooms become visible. 

2.4. Nuclear physics sets the stage 

Nuclear physics therefore governs not only the collapse phase, where electron captures 
on nuclei were recently included in a much improved way and were found to determine 
the position of shock formation and the structure of the layers the shock expands into [ 
132] (also W.R. Hix and G. Martinez-Pinedo, this conference). Nuclear physics and in 
particular the properties of the nuclear EoS also determine the contraction and size of the 
nascent NS and thus may influence the growth of nonradial instabilities and anisotropies 
in the postshock accretion flow during the z^-heating phase. 

In fact, differences in the nuclear EoS have interesting consequences for core collapse, 
bounce conditions, shock formation, and postshock evolution in ID simulations (for 
a brief summary, see [ 121] )• Figure |U shows selected results for three different EoSs 
(Lattimer & Swesty [[Tlj, Shen et al. [EU, and Wolff & Hillebrandt [CSj). The softest 
of them (L&S) leads to the highest densities at bounce and the smallest enclosed mass of 
the shock formation position. Lateron the nascent NS contracts most rapidly, forcing the 
shock to retreat much more quickly than in case of, in particular, the stiff Wolff EoS (Fig. 01 
upper middle panel). These differences affect the v luminosities during the v e burst and 
the postbounce accretion phase (Fig.HJ upper right panel). While for the Wolff EoS the v e 
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Figure 4. Left: Pressure (top) and adiabatic index r = (din P/dln p) s (bottom) vs. mass density 
for an entropy s = 1 k-Q per nucleon (and Y e = 0.4) for the EoSs of Ref. [ 115] ("Wolff", thin lines), 
Ref. [d| ("Shen", medium) and Ref. [HU ("L&S", thick). Middle: Shock positions (solid lines) 
and neutrinospheric radii of v e (dashed) as functions of time for collapse simulations of a 15 Mq 
progenitor (Model sl5a28 [EU) with the three nuclear EoSs (top) and with two different EoSs 
at densities p < 10 11 gem -3 (L&S compared to an ideal gas EoS of e~ , e + , 7's, and Boltzmann 
gases of n, p, a and a representative heavy nucleus in NSE; bottom). Right: Prompt v e burst 
(left panel) and postbounce luminosities of v e (solid lines), v e (dotted) and heavy-lepton z/s and 
P's (dashed) for the simulations of the 15 Mq star with the three different nuclear EoSs (top) 
and the two different low-density EoSs (bottom). (The plots were taken from Ref. [ 33 .) 



release in the burst is highest, a more compact NS and thus hotter neutrinosphere causes 
higher post-burst v luminosities and mean energies. In the lower middle and right panels 
we show the result of a test which we performed with a four-nuclei NSE- EoS replacing the 
low-density (p < 10 n gcm~ 3 ) part of the L&S EoS (where an error in the treatment of 
a-particles has recently been discovered; C. Fryer, J. Lattimer, personal communication). 
The conditions in the postshock layer were hardly affected because the gas is disintegrated 
into free nucleons at the high entropies encountered in simulations with relativistic gravity. 
Minor differences in the evolution of the shock radius (Fig. lower middle panel) and v 
luminosities (Fig. HJ lower right panel) were caused by differences in the mass accretion 
rate associated with the EoS treatment in the infall region ahead of the shock. 
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3. SUMMARY 

Non-radial instability of the accretion shock can amplify vorticity in the postshock 
flow and thus can support the growth of convection in the ^-heated layer. It may be an 
important ingredient for eventually robust explosions by the //-driven mechanism. Low 
modes in the flow can develop when the explosion sets in slowly. The relative importance 
of the different instabilities seems to depend on v cooling and heating on the one hand 
and the high-density EoS, which controls the contraction of the nascent NS, on the other. 
Simulations require the use of a full 180° grid and ultimately may have to be done in 3D. 
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